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We calculate the Coulomb gauge temporal quark and gluon propagators in quenched QCD. From the temporal 
quark and gluon propagators, dispersion relations and quasiparticle masses are determined by means of the 
Maximum Entropy Method. 



It is well known that QCD undergoes a phase 
transition to a deconfined phase where it is be- 
lieved that the dominant degrees of freedom are 
quasiparticles with quantum numbers of quarks 
and gluons. Quasiparticles show up as poles in 
the retarded quark and gluon propagators. The 
corresponding dispersion relations ( the position 
of the poles ) was studied in so-called HTL per- 
turbation theory which is valid if the separation of 
different scales holds 1/T < 1/gT < l/g 2 T (see 
e.g. jl| for a review). In the interesting temper- 
ature region, however, the coupling is large g>l. 
Nevertheless, the corresponding quasiparticle pic- 
ture finds application in refined perturbative cal- 
culationis of the bulk thermodynamic properties 
where it helps to improve the convergence of per- 
turbative series || as well as in more phenomeno- 
logical approaches B. In view of these facts a 
non-perturbative study of quark and gluon dis- 
persion relations in the deconfined phase is highly 
desirable. 

We have calculated the quark and gluon prop- 
agators in quenched QCD with Wilson fermions 
on 64 3 x 16 lattice at two different temperatures, 
T = 1.5T C using 20 configurations and T = 3T C 
using 40 configurations (T c is the critical tem- 
perature of the deconfinement phase transition). 
We have used the standard Wilson action for the 
gauge fields and an 0{a) improved fermion action 
(clover action). Furthermore we have performed 
our simulations at the critical value of the hop- 
ping parameter n = k c (T = 0) determined at zero 
temperature. The values of the gauge coupling 
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corresponding to the two values of the tempera- 
tures are = 6.872 for T = 1.5T C and = 7.457 
for T = 3T C . The non-perturbative value of the 
clover coefficient c sw and of the value of the criti- 
cal hopping parameter k c for these 0- values were 
obtained from nterpolationsof results given in [Q : 
c sw = 1.412488, k c = 0.13495 for T = 1.5T C and 
c sw = 1.338924 , k c = 0.13390 for T = 3T C . 

Lattice calculations can provide information 
on the imaginary-time (Matsubara) propagator 
D{iuj n ,p) (uj n being the Matsubara frequencies). 
This is related to the retarded propagator by an- 
alytic continuation Dr(pq,p) = —D(po + ie,p)- 
This implies 

Dhui n ,p) = - / duj- , (1) 

J-oo ^« - w 

where p(u>,p) = ilm_D^(c<j + ie,p) is the spectral 
function. 

As quark and gluon propagators are gauge de- 
pendent quantities it is necessary to fix a par- 
ticular gauge which we chose to be the Coulomb 
gauge. In this gauge the condition diA{(r, x) = 
is imposed independently on different time slices 
which allows to construct a transfer matrix || 0. 
As the consequence of this the spectral function 
of quarks and gluons is positive which allows to 
use the Maximum Entropy Method (MEM) to re- 
construct the spectral functions (see Ref. for 
a review). Though the quark and gluon propaga- 
tors are gauge dependent quantities the positions 
of the peaks in the spectral function (poles in 

2 To fix the gauge completely an additional time-dependent 
gauge transformation is necessary^ This, however, depends 
on temporal links Uq{t, x) only |6|. 
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Figure 1. The dispersion relation of quarks (a) 
and gluons (b) for temperature T = 1.5T C and 
3T r „ . 



the retarded propagators) is gauge independent 
at any order of perturbation theory ||. Gauge 
independence of the peak position can be proven 
also non-perturbatively in a class of gauges allow- 
ing the construction of a transfer matrix || . At fi- 
nite temperature there are two kinds of quasipar- 
ticle excitations (corresponding to two branches 
in the dispersion relation) Q: the real quasipar- 
ticles (quarks and transverse gluons ), which are 
the analog of partonic degrees of freedom at zero 
temperature, and collective excitations (plasmino 
and longitudinal gluons), which have exponen- 
tially small residues for momenta p > T . 

We have calculated the quark and gluon 
propagators in the mixed (r, ^-representation, 
D(r,p) = J2n D ( iuJ ri,p)exp(-iu} n T) for several 
values of the 3-momenta p. Since the smallest 



non-zero momentum available on our lattices is 
Pmin — 1.57T, the particular behavior of the col- 
lective excitations at small momenta cannot be 
resolved. Using Eq. (|l|) the gluon propagator 
in the mixed representation can immediately be 
written in terms of the spectral function, 



n 9 

(T,L) 



x cosh(w(r - 1/2T)) 
,, sinh(cj/22 ) 



where T and L refer to transverse and longi- 
tudinal gluons respectively. For the longitudi- 
nal propagators at non-zero momentum we have 
found Dl (t, p) = within statistical accuracy of 
our calculations. Using Eq. (|2|) we determine the 
gluon spectral function using MEM and from the 
position of the peak the gluon dispersion relation 
uj g (p). We note, however, that the gluon spectral 
function has a significant continuum contribution 
above the light cone (u> > p) which is absent in 
HTL perturbation theory. 

We also have studied the behavior of the tem- 
poral gluon propagators in terms of the effective 
masses (see e.g. [|| for the definition). Usually 
effective masses reach a plateau which is the zero 
temperature mass or the screening mass at finite 
temperature. This need not be the case for the 
temporal effective masses at finite temperature. 
In fact, the effective masses m g (r,p) extracted 
from the gluon propagators are always larger than 
the position of the peak in the gluon spectral 
function. This can be attributed to the pres- 
ence of the continuum contribution in the spectral 
function. 

The temporal gluon propagators are influ- 
enced by a zero Matsubara mode contribution 
D? T L Aiw n — 0,p), which is the static magnetic 
propagator in momentum space studied in detail 
in Ref. ||. At p — the magnetic propaga- 
tors are strongly volume dependent and very large 
lattices are needed to perform a reliable infinite 
volume extrapolation As a consequence the 
dispersion relation at zero momentum , i.e the 
plasmon frequency uip, cannot be reliably deter- 
mined from our present calculations. 

The most general form of the temporal quark 
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propagator is 

D«(T,p)= l0 F(T,p)+i-nG(T,p) + H(T, P ), (3) 

with ft — p/p. In the chiral limit the last term 
vanishes. Indeed, we have found that H(r,p) = 
within statistical accuracy. Using the most gen- 
eral form for the retarded quark propagator || 
and Eq.(@) one can derive the following represen- 
tation for the functions F and G: 

F(r,p) = du PF ( U ,p) <*^;}gp\ (4) 

G ( T >P) = / °° dLJ PG fa P) ""c^C J/2T) T)) ■ ( 5 ) 

It turns out that our data on G are too noisy 
to apply the MEM analysis to them. The lo- 
cal masses extracted from F and G, however, 
are identical within statistical errors. In fact, if 
the quark propagators are dominated by a sin- 
gle quasiparticle contribution, the local masses 
extracted from F and G should be identical. 
We therefore applied the MEM analysis only for 
F. The reconstructed spectral function pp{uj,p) 
shows only a single peak indicating the absence 
of the plasmino branch for the values of the mo- 
menta studied by us (at zero momentum there is 
no distinction between the quasiparticle (quark) 
and plasmino branch). In contrast to the gluon 
spectral function the quark spectral function ex- 
tracted from F has negligible continuum contri- 
bution above the light cone. As a consequence 
the corresponding local masses m q (r, p) show a 
plateau already at ra = 1 (a denotes the lattice 
spacing). The position of the peak in the spec- 
tral function uj q (p) agrees with the average value 
of the local masses m q (T,p). 

Our numerical results for the dispersion rela- 
tions of quark and gluons uj q (p) and u) g (p) are 
summarized in Fig. 1. While at T = 3T C the 
dispersion relation both for quarks and gluons is 
close to the free dispersion relation oj 2 (p) = p 2 
one sees large deviations from it at T = 1.5T C . 
In order to quantify the deviations from the free 
propagation we have fitted the dispersion rela- 
tions to u]q (p) — p 2 + m q, g an d determined the 
values of quasiparticle masses m q and m g . For 
T = 3T C we have found m q /T = 1.7 ± 0.1 and 
rrig/T — 1.2 ± 0.1. The value of the gluon mass 



is compatible with the leading order result of 
HTL perturbation theory m 2 = g 2 T 2 /2 if we as- 
sume for the gauge coupling g(3T c ) ~ 1.6 as sug- 
gested by the sort distance behavior of the heavy 
quark potential [jioj. With the same value of g 
the leading HTL result for m q is roughly a fac- 
tor 2 smaller than the value found by us. For 
T = l.5T c we have obtained m q /T = 3.9 ± 0.2, 
rrig/T — 3.4 ± 0.3. Here the value of m q was 
obtained by omitting the first three values of 
u>q(p) from the fit. The values of the quasi- 
particle masses at T = 1.5T C are considerably 
larger than the corresponding ones at 3T C and 
those expected from perturbation theory. This 
is consistent with the temperature dependence of 
quasiparticle masses m q>g /T used in quasiparticle 
models for the equation of state || . 
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